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ABSTRACT 



CO and dust emission at millimeter wavelengths are independent tracers of cold interstellar matter, which have seldom been compared 
on the scale of giant molecular clouds (GMCs) in other galaxies. In this study, and for the first time in the Large Magellanic Cloud 
(LMC), we compute the molecular cloud masses from the millimeter emission of the dust and compare them with the masses derived 
from their CO luminosity and virial theorem. We present CO (1 = 1-0 and J=2-l) and 1.2 mm continuum observations of the Nil 
star forming region in the LMC obtained with the SEST telescope and the SIMBA bolometer, respectively. We use the CO data to 
identify individual molecular clouds and measure their physical properties (CO luminosity, size, line width and virial masses). The 
correlations between the properties of the Nil clouds are in agreement with those found in earlier studies in the LMC that sample 
a larger set of clouds and a larger range of cloud masses. For the Nil molecular clouds, we compare the masses estimated from the 
CO luminosity (Xco ^co), the virial theorem (Myii) and the millimeter dust luminosity (Xi.2mra(dust)). The measured ratios Lco/A^vu 
and Xi.2mm(dust)/Mvi, constrain the Xco and 'Xjust (dust emissivity at 1.2 mm per unit gas mass) parameters as a function of the virial 
parameter dvii- The comparison between the different mass estimates yields a Xco-factor of 8.8±3.5 xlO^" cm"^ (K km s"')"' a^j' 
and a'TCdust parameter of 1.5 ±0.5 xlO"^ cm^ g"' ayii. We compare our Nil results with a similar analysis for molecular clouds in the 
Gould's Belt in the solar neighborhood. We do not find in Nl 1 a large discrepancy between the dust millimeter and virial masses as 
reported in earlier studies of molecular clouds in the Small Magellanic Cloud. The ratio between -£i2mm and Myir in Nil is half of 
that measured for Gould's Belt clouds, which can be accounted for by a factor of two lower gas-to-dust mass ratio, as the difference 
in gas metallicities. If the two samples have similar qtvu values, this result implies that their dust far-IR properties are also similar 
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1. Introduction 

Giant Molecular Clouds (GMCs) are the nursery of massive 

stars. They are of prime importance to the evolution of their 

. host galaxies, but there are still uncertainties on how to best es- 

' timate their masses. At their low temperatures (~ 10 - 50 K), 

the molecular masses cannot be measured through the emis- 

■ sion of the H2 molecule because its first rotational transition 

■ arises from an energy level at a temperature of 510 K. Generally 
, in galaxies, low-J transitions from the second most abundant 
' molecule, CO, are the tracer of the molecular mass. Dust emis- 
sion is an alternative tracer of the gas mass, which has been 
used to estimate dust masses dPraine et al.ll2007tlGalametz et alj 
I2OIII) and to calibrate the Xra conversion factor frorn CO lu- 
minosity to H2 mass jlsraej Tl997t iLerov et all 1201 ll) . In the 
Milky Way and nearby galaxies, where GMCs are spatially re- 
solved, the Xco-factor has been also estimated using the rela- 
tionships between cloud properties assum ing virial equilibrium 
dSolomon et alJll987tlBolatto et al.ll2008l) . 

CO emission does not have a one-to -one correspo ndence 
with H2. Observations in the Galaxv (|& enier et al. 120051: 
lAbdo et al.ll2010l: iPlanck Co llaboration et al."201 Id) and mod- 
els of photodissociation regions (PDRs) (e.g. Wolfire et al. 2010) 
show that CO does not trace the molecular gas over the full 
extent of clouds. CO-dark H2 gas was also observed in the 



Large Magellanic Cloud (LMC) with dus t observations from 
Spitzer a nd Herschel jBernard et al. 2008; iRoman-Duval et all 
2010; GalHanoetal. 2011). In the Small MagellanicCloud 
(SMC). .Rubio et al.. (2004) and later iBot et al.i (iJOOTllIoToh re- 
ported differences up to a factor ~10 between GMCs masses de- 
rived from the dust emission (based on 870 fim and 1.2 mm 
emission) and the virial theorem. Additional evidence for H2 
gas without CO emission is provided by observations of the 
[Cii]/ll58//m l ine em ission in Magellanic ii Tegular galaxies 
(Madden et a l] 119971; [Israel & MalonevI 1201 lb and the Milky 
Way (iVelusamv et alJl20 10). PDR models show that the column 
density threshold for CO to become the main carbon species has 
a much stronger dependence o n the metallicity tha n that for the 
gas to become molecular (e.g. I Wolfire et al.ll2010h . This is be- 
cause CO is shielded from photodissociation by dust while H2 
is self-shielded. Whereas in the Galaxy ther e is a c ontribution 
to CO emission from diffuse molecular gas (iLiszt et al. 201 (J), 
in low metallicity galaxies the CO may predominately come 
from dense, high column density clumps embedde d in a diffuse 
molecular envelop w ith weak or no CO emission ([Israelii 19881 
lLequeuxet"anil994 . 

The motivation of this paper is to compare estimates of the 
molecular cloud masses in a l ow metallicity envir onment, the 
LMC. The LMC is at 50 kpc dPersson et ai1l2004 . It presents 
a metallicity, measured for Hit regions, of 12-i-log(0/H)=8.37 
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(iKeller & Woodl2006t) . that corresponds to Z^O.5 Zq. Extensive 
CO(l-O) surveys have been done in the Magellanic Clouds. In 
the LMC, the first one done by Cohen et al. (1988) a nd sub- 
seque ntly the SEST (Israel et al. 1993), NANTEN (Fuk uietal.1 
l2008h and MAGMA (.Hughes et aL,2010:,Wong et aL,20lT l) sur- 
veys gave a large set of GMCs. Also, several deta iled studies 
have been done in the LMC (e.g. Garav et al. 2002 for molecu- 
lar complex No. 37; Israel et al. 2003a for Ni l; Johansson e t al] 
[19981 iKutner et all [19971 iPineda et al J 120091 for 30Dor). All of 
these studies report a CO-to-Ha conversion factor estimated from 
the comparison between the virial and CO luminosity masses of 
the molecular clouds larger than in the Galaxy. 

This paper is based on CO(J=1-0 and J=2-l) and 1.2 mm 
continuum observations of the Nil star- forming region in the 
LMC. We present the observations and compare the CO lumi- 
nosity of molecular clouds with their masses estimated with the 
viiial theorem and the millimeter dust emission. Unlike what has 



it is observed that generally [C n] does not peak where CO peaks 
but it is offset in the direction of the ionizing stars. 



been done with far-IR obser vations from Spitzer and what is now 
being done with Hersch el ([Bernard et al.l 120081; JGalhano et alj 
1201 it JLerov et aDl201 i h. the computation of molecular masses 
from the millimeter emission of the dust is not very sensitive to 
the dust temperature because it is on the Rayleigh- Jeans spectral 
side of the energy distribution. Further, the comparison between 
dust and virial masses has seldom been performed because the 
molecular clouds have to be resolved to estimate virial masses. 
With single dish telescopes and space IR observations (Spitzer 
and Herschel) this is possible only in the Milky Way and the 
Magellanic Cl ouds. This p aper extend s the earli er studies on the 
SMC done bv lRubio et al] (2004) and JBotetal . (2007, 2010) to 
intermediate metallicities in the LMC. Soon with ALMA, this 
can be extended to more distant galaxies in the Local group. 

The paper is organized as follows. In Section 2, we present 
observations of CO in the J=2-l transition line, new data in the 
transition J=l-0, and a continuum map at 1.2 mm. We identify 
the Nil molecular clouds using the CO(2-l) data and derive 
their physical properties (CO luminosity, size, line-width and 
virial mass) in Section 3. In Section 4, we derive millimeter dust 
fluxes from the continuum observations. In Section 5, we present 
the three estimations of molecular masses, and the comparison 
between them. We show for which values of the Xco-factor and 
the millimeter dust emissivity per hydrogen, as a function of the 
virial parameter of the clouds, the masses are identical. Finally, 
Section 6 presents the conclusions. 



2. Observations in tlie N11 region 

2. t. The N11 star-forming region 



Nil (iHenizd [ 19561) is the second largest and brightest Hii 
region in the LMC, and hosts several stellar clusters (see 
Fig. [Tj). Its main feature is a cavity with an inner diameter 
of ~7 00" (~170 pc). In its center there is the OB association 
LH9 (L ucke & Hodgd [19701) while in its northeast part there 
is the NllB nebula which is associated with another stellar 
cluster, LHIO. The c l usters are 7 and 3 Myr old, respectively 
("Mokiem et al."2007h. [Parker et al.l (Il992l) show that most of the 
ionizing flux in N il come s from these two clusters. IR obser- 
vations (Barba et a l. 2003) r eveal th at the star formation is on- 
going in NllB. [fsraeletal.l (l2003ah presented CO(l-O) obser- 
vations of Nil with a resolution of FWHM=45". They identi- 
fied molecular clouds and computed their physical properties. 
Recently, Israel & Maloney (20 11) showed that even though the 
[Cn]/il58//m and CO emission lines present similar spatial dis- 
tributions, the [C n] emission extends beyond that of CO. Also, 
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Fig. 1. Image obtained from the combination of hydro- 
gen, sulfur and oxygen bandpasses (Magellanic Clouds 
Emission-Line Survey MCELS, C. Aguilera, C. Smith and S. 
Points/NOAO/AURA/NSF). Contours correspond to the CO(2- 
1) line emission. We indicate the position of the OB associations 
in Nil, LH9 in the central cavity of Nil, LHIO in the nebula 
NllB, LH13 in Nl IC and LH14 in Nl IE. 



2.2. SEST Observations: CO emission lines 

The observations of the CO emission lines were taken with 
the SEST (Swedish-ESO Submillimeter Telescope) at La Silla 
Observatory, Chile, during January and September in 2001. 
SEST was a single dish antenna of 15 meter of diameter op- 
erating in the range of frequencies 70-365 GHz. It was decom- 
missioned in 2003 with the beginning of operations at the APEX 
(Atacama Pathfinder Experiment) telescope. 

We mapped the central part of the Nil star- forming region 
in three different regions, which are listed in Table [T] We ob- 
served simultaneously the CO J=l-0 and CO J=2-l lines at the 
frequencies 115 GHz and 230 GHz with a full width beam at 
half maximum (FWHM) of 45"(1 1 pc linear size) and 23"(6 pc 
linear size), respectively. The spatial sampling of the telescope 
pointings for both frequencies is 24". We note that our CO(2-l) 
observations are not fully-sampled. However, this does not have 
a significant impact for the goal of this paper (see § 13.1) . The 
system temperature varied between 360 and 713 K at 115 GHz 
and between 235 and 467 K at 230 GHz, during the observations. 

The observations were done in position switch mode with 
a fixed reference off' position free of CO emission. A narrow 
band AOS high-resolution (HRS) spectrometer with 2000 chan- 
nels, 80 MHz bandwidth and 41.7 kHz channel separation (cor- 
responding to 0.055 km s ' for the '^CO J= 2 - 1 line) was used 
as back end. Intensity calibration was done using the standard 
chopper-wheel technique. The pointing accuracy, checked dur- 
ing the observations on RDor, was better than 2". 
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Fig. 2. Integrated CO(2-l) emission line of the Nil star forming region, observed with the SEST. Contour corresponds to 2.5 
K km s~' . The dots represent the observed positions. Red crosses mark the position of the molecular clouds listed in Table|2] Offset 
positions are relative to a :04''56'"42'.47, 6 :-66°28'42'.'0 J2000.0. 



Table 1. Characteristics of the SEST observations in the Nil 
star-forming region. 



Zone 



Coordinates 
R.A. DEC. 



K 



Map! 04"57'"25' 
Map 2 04''55'"36' 
Map 3 04''56"'25' 



.00 
.00 
.00 



-66°24'0a'0 
-66°30'0a'0 
-66°36'0a'0 



0.25 
0.22 
0.20 



0.30 
0.33 
0.27 



" computed within the 0.055 km s ' channel width. 
' computed within the 0. 11 km s"' channel width. 



The data reduction was done with CLASS (Continuum and 
Line Analysis Single-dish SoftwareQ). The complete set of spec- 
tra were visually inspected one-by-one, we did not find any 
anomalies in the spectra. We fitted and subtracted a baseline 
of second order to each spectrum to correct for residual at- 
mospheric and instrumental emission. The observed velocity 
resolution is 0.11 and 0.055 km s-\ for CO(l-O) and CO(2- 
1), respectively. Finally, to obtain the radiation temperature of 
the sources, we corrected for the antenna effic iency ri-Q.l and 
0.5 f or CO(l-O) and CO(2-l), respectively ([Johansson et al.l 
Il998h . The flux measurements have a photometric calibration 
uncertainty of 20%. Table [T] lists the Icr noise in a 0.11 and 
0.055 km s ' channel band, for CO(l-O) and CO(2-l) respec- 
tively, in each line emission and mapped zone. These values 



http;//www.iram.fr/IRAMFR/GILDAS 



were estimated as the average value of the standard devia- 
tion computed in 20 line-free channels at velocities of about 
259kms-i. 

Figure |2] shows the integrated CO(2-l) line emission Iqo, in 
units of K km s"'. This image was obtained by integrating the 
CO emission in a velocity window defined in the total spectrum, 
from 265 to 289 km s ' . Figure |3]pi"esents the CO(2-l) channel 
maps integrated in bins of 2 km s"'. From these channel maps 
we can distinguish several components at different velocities. 

For each CO transition line, we computed the total emission 
by adding all the spectra in the field and integrating the emission 
in the velocity range 265-289 km s"' . The total emission corre- 
sponds to 1416 K km s' and 1630 K km s' for CO(l-O) and 
CO(2-l), respectively. 

2.3. SIMBA observations: millimeter continuum emission at 
1.2 mm 

To measure the dust emission in Nil, we observed the con- 
tinuum emission at 1.2 mm. This observation was obtained 
with the SEST IMaging Bolometer Array (SIMBA) at La Silla 
Observatory in October 2002. SIMBA was a 37 channel hexago- 
nal array operating at 250 GHz. The bandwidth in each channel 
was about90 GHz. The Half Power Beamwidth (HPBW) of each 
element was 24". 

The observations covered the central and northern parts of 
the Nil complex. The individual maps were produced using 
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Fig. 3. Channel maps for the Nil region. Contours illustrate the CO emission in the transition J=2-l smoothed over 4 pixels. The 
first contour is 0.75 K km s"' , with a spacing of 0.8 K km s"' . From the top-left to the bottom-right image, we present channel maps 
starting from 266 km s"' to 290 km s ', with a velocity width of 2 km s '. Offset positions are the same as in Fig. 12.21 We mark in 
red the molecular clouds identified with CPROPS. 



the fast scanning mode. The elimination of the correlated sky 
noise, the co-addition of individual maps, and the photometry 
was done using the MOPSI packagqj The t zenith opacity was 



^ MOPSI is a data reduction software package for IR and radio data 
developed by R. Zylka, IRAM, Grenoble, France 



determined performing sky dips every 2 hours during the day, and 
every 3 hours during the night. The flux determination was done 
for each individual observing run by observing Uranus. The flux 
measurements in the final coadded map have a photometric cal- 
ibration uncertainty of 20%. The sensitivity of the final image 
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is 4 mJy beam ' (Icr). Figure |4] shows the observed emission of 
Nl 1 at 1.2 mm, with CO(2-l) contours overlaid. 



' 'J J " -A 




-^^ 





04'^ 58' 



■ if-- ■ ■■ I ■ ■ I 



152 



128 



103 



79 



54 



30 



U5 



57 56 

Right Ascension 



Fig. 4. Millimeter continuum emission at 1.2mmof theNll star 
forming region, observed with the SIMBA. Black contour rep- 
resents the emission at 3o--12 mJy beam"'. Red contour is the 
CO(2-l) emission at 2.5 K km s"'. 

This image was used to measure the millimeter continuum 
flux of the molecular clouds in Nil, identified with the CO ob- 
servations. We measure the emission in an elliptical aperture de- 
fined by the molecular cloud extent at the position of the central 
coordinates of the cloud. We subtracted the nearby sky emission 
measured next to the source on an empty image area. The flux 
errors correspond to the noise level (standard deviation) in units 
of Jy beam ', measured locally on an OFF area next to the cloud, 
multiplied by the area of the beam and by the root square of the 
number of beams within the aperture size. The final flux errors 
were obtained using error propagation of the flux uncertainties 
and the 20% photometric calibration errors. 

2.4. Ancillary data 

We complement our data set with two published fully reduced 
images at 8.6 and 4.8 GHz obtained with the ATC A (Australia 
Telescope Compact AiTay) in February 2003 (Dickel et al. 
[2005). These observations are part of two surveys of the LMC. 
The angular resolution is 22"and 35"(FWHM beam width) at 
8.6 and 4.8 GHz, respectively. The flux uncertainties for the ob- 
servations are 5%. To measure the fluxes of these radio emission 
images, we did aperture photometry in the same way as for the 
emission at 1.2 mm. 

3. Molecular clouds in N11 

In this section, we present the observational results obtained 
from the SEST observations in Nl 1. 

3.1. Identification of molecular clouds 

The identification of the molecular clouds in Nl 1 was done us- 
ing the CO(2-l) line emission observations, which has an angu- 



lar resolution twice that of the CO(l-O) line emissi on. We use 
file CPROPS algorithm (iRosolowskv & Lerovll2006h to identify 
the individual molecular clouds. This program has been used in 
sever al other studies of molecular clouds, i ncluding the LMC 
(e.g. iFukui et al.ll200i iHughes eTaDIMol: IWong et al.ll201 H) . 
CPROPS is an IDL procedure which identifies clouds from a 
data-cube using the spectral and spatial structure of the line 
emission. The physical parameters of the identified molecular 
clouds are measured as moments of the emission along the spec- 
tral and spatial axis. The CPROPS algorithm works as follows. 
It defines a mask from a cr,-ms value of the data. Within this mask, 
pixels with emission higher than a threshold value of 5cr„„s, over 
at least two consecutive velocity channels, are identified as cores 
of the mask. The mask includes pixels with brightness larger 
than IcTrma emission which meet the two following conditions. 
They are spatially and spectrally connected to the core pixels 
and they are not spectrally isolated (i.e. at least two consecu- 
tive channels are above 2crr„s). Every core is decomposed into 
candidates molecular clouds. The parameters used for this de- 
composition were: minimum area of 4 beams, line-widths larger 
than 1.5 km s"', and a minimum ratio between the flux of the 
cloud at the peak and at its edge of 3o-^ms- The algorithm defines 
elliptical molecular clouds, which are extrapolated to the K 
isosurface to estimate what we would measure with no noise. 
For a detailed explanation of CPROPS, we direct the reader to 
Rosolowskv & Lerov (2006). 

To increase the SNR in our data, we smoothed the data in 
the velocity axis to a spectral resolution of 0.5 km s ' . CPROPS 
identified 22 individual molecular clouds. We discarded three 
clouds with low signal-to-noise ratio (<4) and which areas are 
close to our 4 beams threshold. Our CO(2-l) data is not fully- 
sampled. It was observed at about twice the Nyquist frequency 
(§ I2.2| |. However, the undersampling of the observations do not 
impact the conclusions of this paper. We quantify the impact of 
the undersampling on the parameters of the clouds identified by 
CPROPS with a simple test. We run CPROPS for our fuUy sam- 
pled CO(l-O) (angular spacing of 24" and beam size of 45") 
data cube with a minimum cloud area of one CO( 1-0) beam. The 
results of this run can be directly compared with the CPROPS 
run on the CO(2-l), for which we have used a minimum cloud 
area of four CO(2-l) beams. In Figure |5] the cloud radii and 
velocity dispersions for both runs on the CO(l-O) and CO(2-l) 
data-sets are compared. The eiTorbars are 1-sigma eiTors on the 
parameters computed by CPROPS. Although the two cloud pa- 
rameters do not match exactly one-to-one, the two sets of points 
are statistically close to each other. There is no systematic diff'er- 
ence in radii, and only a slight difference by a factor of 1.2 be- 
tween the mean velocity dispersion for the CO(l-O) and CO(2- 
1) lines. We conclude that the undersampling of the CO(2-l) 
observations has a minor impact on the paper conclusions which 
are all of statistical nature. The difference in line widths implies 
a factor 1 .5 systematic on our estimates of the virial masses. We 
do not correct our 2-1 line widths for this factor but consider it 
as a systematic uncertainty in our scientific analysis. 

Figure |2] presents the positions of the molecular clouds in 
the Nil region. Channel maps in Figure |3] show the boundaries 
of the clouds. Table |2] lists the final set of individual molecular 
clouds identified with CPROPS and their properties. The radii 
correspond to the geometrical mean of the ellipse semi-axis. 
CPROPS gives velocity dispersions which we convert to line- 
widths at half maximum with the relation Av = cr^, 2 y/l ln(2). 
Uncertainties on the radii, line-widths and CO luminosities come 
from the fitting procedure in CPROPS (IRosolowskv & Lerovl 
2006i) . Error bars on the CO luminosity also include the 20% un- 
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Fig. 5. Comparison between the radius and velocity dispersion, 
for molecular clouds obtained from CPROPS runs the CO(l-O) 
(red squares) and CO(2-l) (black triangles) data set. The error 
bars correspond to the 1-sigma errors computed by CPROPS. 



certainty on the calibration. Clouds #1 and #5 present different 
velocity components and might be separated into 2 clouds each. 
The mean local standard of rest velocity Vlsr is 278 km s"'. 

The total intensity of the identified molecular clouds corre- 
sponds to 61% and 50% of the total emission (see § I2.2| | for the 
CO(l-O) and CO(2-l) emission lines, respectively. 



Ref. 

# 



This work 



Israel et al. (2003a) 
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2.8 
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10 


13 
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15 


14 


10.4 
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11.2 
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15 
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4.2 


25.2 
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11.3 
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18 
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26.2 
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17 
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24 
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Table 3. Comparison between molecular clouds detected 
in CO(2-l) and previous clouds detected in CO(l-O) by 
Israel et al. (2003a ). To compare the masses, we scale the Myir 
values from Israel et al.l (12003 a) to use the formula for spherical 
clouds with a density profile oc r~' . 



3.2. Comparison with previous data sets 

The Nil star- forming region was previously observed in CO(l- 
0) as part of the ESO-SEST Key Programme (Israel et al. 2003a) 
with an angular resolution of 45". They identified 14 molecular 
clouds in the same area that we cover. We compare our list with 
their results. All of the clouds previously identified are in our 
list of clouds. Our CPROPS run could not resolve two of their 
clouds (#16 and #17, Tabled, due to the undersampling of the 
observations. These clouds are small (~30", i.e. 7 pc, of radius) 
and close to each other (40", i.e. 10 pc, less than twice the sam- 
pling of our observations). Thanks to the higher J=2-l CO reso- 
lution, we find six additional clouds. Table|3]lists the correspon- 
dence between both studies, comparing radii, line-widths and 
virial masses. Even though there are differences in the data-sets 
and methodologies, global parameters of the molecular clouds 
ar e similar. St a tistical ly, our line-widths are larger than those 
of Israel et al. (2003a), which yields larger viiial masses. We 
checked with the data that this difference mainly co mes from the 
way the line-widths are measured. Those reported by llsrael et al] 
(l2003ah are measured on the spectrum at the emission peak of 
the molecular clouds, while our values are computed on the spec- 
trum integrated over the entire molecular clouds. The higher sen- 
sitivity and resolution of our data could also contribute to this 
difference. We may have included weaker velocity components 
located in the same line-of-sight of the clouds, which were not 
detected before. 

We also compare our results with those of the MAGMA sur- 
yey in the w hole LMC (CO(l-O) observations, FWHM=45") 
dWong et al.l l201 lb. Like us, they identified molecular clouds us- 
ing CPROPS. The parameters used in their decomposition to 
identify molecular clouds are a minimum ratio between the flux 
of the cloud at the peak and at its edge of 2cr. In the same area 
that we observed, they identified 18 clouds. 15 clouds in their 



and our list coincide. Two of the remaining three clouds are part 
of a single cloud in our list, while the third one is small (<5 pc) 
and was not detected by our CPROPS run. Statistically, our 
CO(2-l) molecular clouds present line-width slightly broader, 
by a median factor of 1.3. If we take into account the 1.2 factor 
which comes from the undersampling of the CO(2-l) observa- 
tions - offset value between the mean velocity dispersions of the 
CO(l-O) and CO(2-l) clouds presented in Section [ll] -, this 
factor increases slightly to 1.5. A comparison of our molecular 
clouds with those identified across the entire LMC ( Wong et al.l 
1201 Ih is displayed in Figures |7] and |8] which show that the Nil 
star-forming region is a singular region in the LMC. 

3.3. Line ratio R2-1/1-0 = /co(2-i)//co(i-{)) 

As described in §2, we mapped the central part of Nil simul- 
taneously in CO(l-O) and CO(2-l) line emission. From these 
lines, we compute the ratio of line intensities in K km s"', 
/?2-i/i-o- To compare these lines, we need to have the data at 
the same spatial resolution. Thus, we convolve the CO(2-l) data 
set (FWHM=23") with a gaussian kernel to the 45" resolution of 
the CO(l-O) data. Maps of integrated emission for CO(l-O) and 
CO(2-l) are obtained similarly as described in § 12.2] by integrat- 
ing the emission within a spectral window in the velocity range 
265-289 km s"^ From these two maps, we derive a /?2-i/i-o ratio 
map for the entire region. We estimate the average /?2-i/i-o ratio 
for each molecular cloud by measuring the emission within an 
aperture defined by the size of the molecular clouds as identified 
by CPROPS. Uncertainties on this ratio include the 20% photo- 
metric uncertainties. Table |2] lists these values. Figure |6] shows 
the CO(l-O) and CO(2-l) intensity relation for each pixel in our 
observations of the Nil region and for each molecular cloud. 
The median value for the /?2-i/i-o ratio is 1.0. This is close to 
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Table 2. Molecular clouds in Nil identified in CO(2-l) with CPROPS. Offset positions are relative to a :04''56"'42".47, 
6 :-66°28'42'.'0 J2000.0. The radii correspond to the geometrical radii for each elliptical cloud. Uncertainties include the error- 
bars in the CPROPS fitting and 20% error in the luminosity due to calibration. "Clouds unresolved by CPROPS, we use the beam 
size as upper limit. 'Molecular masses estimated from the CO(2-l) luminosity, scaled to CO(l-O) by the /?2-i/i-o values, using 
the Zco-factor for the LMC estimated by Hu ghes et al.1 (1201 Ol) of 4.7 X 10^°cm"^ (K km s~' r' (see section ISTTTl. '| Virial masses 
computed with the formula Mvii = 190 Av^ R for a spherical cloud with a density profile oc r'^ dMacLaren et al.|[l988l) . 



the median value of 1.2 reported by llsrael et al.l (l2003bb for an- 
other sample of molecular clouds in the LMC. In Figure |6] we 
observe some dispersion of the pixels around the line. There are 
a few outliers. Four molecular clouds, 2, 6, 8 and 10, show large 
^2-1/1-0 values. Clouds 6 and 8 are associ ated with th e embed- 
ded IR stellar cluster in Nl 1 identified bv iBarba et al.l (12003.) . 
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Fig. 6. /co(i-O) vs. /co(2-i) line emission for each pixel in Nil 
(black diamonds) and for each molecular cloud in Table |2] (red 
circles). The cloud fluxes were divided by 7. The line corre- 
sponds to the median value of the ratio /?2-i/i-o=l-0. 



3.4. Correlations between physical parameters 

In this section we compare the physical properties of the Nil 
clouds with empirical relations derived from wider surveys in 
the Galaxy and the LMC. We determine the correlations between 
the physical properties of the molecular clouds Lco(i-O), R, A^vir 
and (Xv presented in Table |2] Lco(i-O) values are estimated as 
the Lco(2-i) measures scaled by the ^2-1/1-0 estimates. Since the 
range of the molecular clouds luminosities in Nl 1 is limited, we 
cannot estimate reliable values for the exponents of the power 
law correlations between the different paramete rs. Thus, we fit 
the N l 1 relations with the exponents found by ISolomon et al] 
(Il987h for GMCs in the Galactic molecular ring. The multi- 
pUcative factor of the following equations result from an error- 
weighted fit of all data points and the error bars are the dispersion 
of the individual points. 



Leo = 6.9^1) R 



+9.0 n2.5 



Leo - 42_28 



^^' crl 



- Q1+79 r0.81 

- ^^-42 ^CO ' 

cr, = 0.72t;^/?«^ 



(1) 
(2) 
(3) 
(4) 



Figure Q shows the relations Leo vs. cr,,, Leo vs. R and o",, 
vs. R. Figure |8]presents the relation between the virial mass and 
the CO luminosity. In these figures, we compare the Nil molec- 
ular clouds with LMC clouds identified by Wong et al. (20lll) 
as part of the MAGMA survey (green crosses). The best fits 
for our clouds, using the Galactic exponents, correspond to the 
solid lines and the shadow areas to the Icr-error in the multi- 
plicative factor. We also include, in segmented lines, the rela- 
tions obtained by Solomon et al. ( 1987) for the Galactic molec- 
ular clouds, Bolatto et al. ( 2008) for GMCs from dwarf galaxies 
and lHughes et aL (.2010 ) for GMCs in the LMC. 
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Fig. 7. Correlation graphics, a) log(Lco) vs. o-y, b) log(Lco) vs. R 
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for individual data points. Black dot-dashed lines: Galactic re- 
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dwarf galaxies from iBolatto et al.l (l2008l) . Green crosses repre- 
sent the molecular cloud s iden tified in the MAGMA survey of 
the LMC by Wong et al] (l20Tlh . 



4. Millimeter dust emission at 1.2 mm 

In this section, we estimate the emission of the dust in 
Nil at 1.2 mm combining the SIMBA (FWHM=24"), SEST 
(FWHM=23") and ATCA 8.6 GHz (FWHM=22") observations. 
To do that, we convolved the SEST and ATCA images to the spa- 
tial resolution of the SIMBA observations. 

In a molecular cloud, the continuum emission at 1.2 mm. 



, is the sum of free-free and dust radiation. Our observa- 



o total 

'^ 1.2mm 

tions obtained with the SIMBA bolometer (see § 12. 3b provide us 
with the continuum emission at 1.2 mm. The CO(2-l) emission 
line falls in the SIMBA bandwidth. Thus, to measure the dust 
emission at 1.2 mm for the Nl 1 molecular clouds, the free-free 
and CO line contribution have to be subtracted. 



We use centimeter-wave observations (at 8.6 GHz) to esti- 
mate the contribution from free-free to the 1.2 mm emission. 
Continuum emission at centimeter wavelengths is composed 
by two contributions, thermal emission from H ii gas (free-free 
emission), and synchrotron emission. Since the spectral index of 
the thermal emission is flatter than that of the synchrotron emis- 
sion and the synchrotron emission is not correlated with GMCs, 
we assume that the emission at 8.6 GHz is mainly tracing the 
free free emission. We scale the radio continuum emission at 
8.6 GHz to 250 GHz (1.2 mm). For an optically thin medium, as 
the emission at sub-mm/mm wavelengths, the free-free emission 
depends on the frequency as 5v oc v""'. From this relation, we 
estimate the free-free emission at 1.2 mm as: 



nff 

'^ 1.2mm 



250 GHz y 
8.64 GHz I 



-0.1 



8.64 GHz- 



(5) 



Uncertainties on S'j^jmm ^^ computed as the quadratic sum of 
the errors, including the flux calibration uncertainty (5% of the 
flux in the convolved and scaled image), and the uncertainty on 
the free-free contribution. The latter was estimated comparing, 
for each cloud, the emission at 4.8 and 8.6 GHz smoothed to 
a common resolution of 35". The error is the absolute value of 
the difference between the 1 .2 mm free-free flux estimated from 
the 8.6 GHz observation and that computed from the geometric 
mean value of the 4.8 and 8.6 GHz fluxes. 

We estimate the contribution of the CO(2-l) emission line 
to the SIMBA bolometer as follows. This contribution can be 
written as 5'co(2-i)=2kv-'c"-'f2/co/Av where v=230 GHz is the 
frequency of the CO(2-l) line, Ico the emission line intensity, Q. 
the solid angle of the Gaussian beam of SIMBA in steradian, and 
Av=90 GHz the frequency width of the bolometer. We obtain an 
expression that we use to compute the contribution of the CO(2- 
1) line and its error bar: 



cCO 

'^ 1.2mm 



mJy beam 



0.20 



tco 



Kkm s 



(6) 



To obtain the dust emission, we subtract the free-free and CO 
contributions from the total continuum emission at 1 .2 mm. The 
error on the dust emission is the quadratic sum of three errors: 
the error on the 1 .2 mm flux including the 20% uncertainty on 
the photometric calibration, the error on the free-free subtrac- 
tion and the error on the CO line contribution. Figure |9] shows 
the different contributions to the emission at 1.2 mm in the Nil 
region. From the left to the right images: the total emission at 
1.2 mm, the free-free emission at 1.2 mm estimated from the 
ATCA emission at 8.6 GHz, the CO(2-l) line contribution to 
the total 1 .2 mm emission observed with SIMBA, and the result- 
ing dust emission. Table |4] lists the fluxes contributions to the 
total 1.2 mm emission for each molecular cloud, measured in 
the images (convolved to the same spatial resolution) by aper- 
ture photometry using the molecular cloud sizes computed by 
CPROPS and subtracting the nearby background within an aper- 
ture of similar size. 



5. Molecular masses 

In this section, we discuss and compare three methods to esti- 
mate the masses of the molecular clouds based on the CO lumi- 
nosity, the virial theorem and the millimeter dust emission. 
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Fig. 9. Different contributions to the emission at 1.2 mm. a) Total emission observed with SIMBA. b) Free-free emission at 1.2 mm 
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images are in units of mJy beam '. Intensities of the images b) and c) were amplified by a factor 6. In all the images, the black 
contour correspond to 3cr (12 mJy beam"') of the total 1.2 mm emission. 



5.1. Masses from CO luminosity, Mqo 



In the Galaxy, the mass of a molecular cloud is tradition- 
ally estimated from its CO emission. The ratio between the 
CO(l-O) integrated intensity Wco and the H2 column density 
A^(H2), is the Xco-factor (A^(H2) = Xco Wco)- This conver- 
sion factor is uncertain because it depends on the cloud prop- 
erties. In the MiUcy Way, it has b een calibrated with y-ray 
emission (i.e. iBloemen et alj 19861). the assumption of virial 
equilibrium (i.e. Solo mon et al.|[l987l) . and dust extinction (i.e. 
^oinbardi et al. 2006). These empirical studies yield values in 



the range 2-4x10 cm (K km s ) . Most recent v alues from 
y-ray in the Gould Belt are lower (lAbdo et al.ll2010l) . but in av- 



erage in the Milky Way the studies agree with the previous range 
of values (Bolatto et al. 2013). 



Several studie s (e.g.. ICohenetal.1 Il988t [Israeli 119971; 
iFukui et al.1 120081: iHughesetalJ l2010h have shown that the 
Galactic conversion factor does not apply to molecular clouds 
in the LMC, where the estimated values are between 2 and 4 
times higher than the Galactic one. In this work, we use the 
Xco-factor computed from the molecular clouds observed and 
identified by the MAGMA survey in the LMC JHughes et alj 



2010) which con-esponds to 4.7 x lO^^cm"^ (Kkms"')" ■ To 
estimate the masses based on the CO luminosity (Mco) by us- 
ing the Xco(i-O) factor estimated from the MAGMA survey, we 
scale the Lco(2-i) luminosities by the /?2-i/i-o values, both listed 
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in Table |2] The inferred CO luminosity masses are also listed in 
Table|2]and they range from 2 to 80 xlO^ Mq. 

5.2. Masses from the virial theorem, My^ 

The virial masses are computed from the virial equilibrium be- 
tween gravitational and turbulent kinetic energies. This implies 
that the molecular clouds are gravitationally bound, ignoring ex- 
ternal pressure and magnetic energy. For a spherical cloud with a 
profile density of p oc r"', the virial mass is Myir - 9/2 Rcrl/G, 
where R is the radius, cTi, the velocity dispersion of the gas in 
1 di mension, and G the gravitational constant ( JMacLaren et aLJ 
[T98&) : 



[Mq] 



190 



Av 



R 



[pc]/\[kms"'] 



(7) 



To compute Mvir we use the line-width and the size estimated 
from the CO(2-l) emission, both listed in Table|2] The resulting 
Mvii- values are listed in the same table. These masses vary be- 
tween 7 and 130 xlO^ Mq. The Icr uncertainty corresponds to 
the quadratic sum of the errors on the line-width and the radius. 
In Table |2l we also list the ratio between Mvii and Mqq. The me- 
dian value for this ratio is 1.9 with a median absolute deviation 
of 0.9. This ratio increases to 2.8 if we include the offset in the 
velocity dispersion due to the undersampling (§ 13.1) . The me- 
dian value of the ratio between the virial mass and Lco(i-O) is 21 
Mq (K km s"' pc^)"' (31 Mq (K km s"' pc^)"' including the ef- 
fect of the undersampling on the velocity dispersion). The ratios 
for all clouds, but clouds #14, 18 and 19, are within 3cr of this 
median ratio. 

5.3. Masses from the millimeter dust emission, M^m 

The masses of the molecular clouds can also be estimated from 
the millimeter emission of the dust. In this method, we assume 
that the dust properties in Nil are the same as those in the 
Galaxy. There is still an open debate about this hypothesis, since 
the fa r-IR dust properties i n the Gal axy and in the LMC may 
differ dGalhano et al.ll201 UlPlanck Co llaboration et al. 201 Ic). 

We approximate the spectral dependence of the dust emis- 
sion with a grey body. Assuming that the emission is optically 
thin, the dust optical depth t^(v) can be written in terms of the 
dust column density. Then, the gas mass estimated from the dust 
millimeter emission is; 



SyD^ 






(8) 



with S V the millimeter cold dust emission, D the distance to the 
LMC, By the black body emission at a dust temperature T^, k^ the 
absorption coefficient per unit dust mass and xj the dust-to-gas 
mass ratio. We can estimate k^ from the emissivity per hydrogen 
euiv) as Kdiv)-eHiv)(xdfimuy^ , where fi is the mass per hydro- 
gen atom of 1 .36. Our masses derived from dust emission, like 
the masses derived from CO luminosity (in section ISTt and the 
virial masses (in section |S!2] |. include the contribution from he- 
lium ji = 1.36. We estimate eH(l-2mm) using the recent results 
obtained with the Planck satellite. By combining Planck and 
IRAS o bservations in the Galaxy, iPlanck Collaboration et alj 
(1201 Ibh fitted the spectral energy distribution of the local ISM 
with a grey body function, finding rciust=17.9 K and j6=1.8. 
They derived a dust emissivity at 250 ^ m which agrees wit h 
that obtained from COBE observations dBoulanger et al.lfl996h . 



From these Planck results, we compute the dust emissivity 
per hydrogen atom at 1.2 mm as eH(l-2mm)=6.1xlO"^^ cm^. 
Furthermore, Planck Collaboration et al. (2011a) showed that 
the emissivity T25()^,„/A^h increases by a factor of two from the 
atomic to the molecular gas in the Taurus molecular cloud. 
Thus, for the dust-to-gas mass ra tio in the solar neighborhood 
of Xd = 0.007 (iDraine & Lill2007h . the absorption coefficient for 
dust in molecular clouds i s /Cd(1.2mm)=0.7 7 cm^ g"' . This value 
is similar to that used by iBot et al.) (l2007h to compute the mil- 
limeter masses of the molecular clouds in the SMC. 

To compute the millimeter masses, we use a single value 
of 20 K for the dust temper atures. This mean value is derived 
from the temperature map of iPlanck Collaboration et aP (1201 id 
their Figure 7) at the position of the Nil region. The same figure 
shows that, around Nl 1, the dust temperature ranges between 15 
and 25 K. We will use this range as the maximum uncertainty on 
the temperatures. 

We assu me that th e dust- to -gas mass r at io sca les with 
metallicity dlssa et alj 119901: iJames et alj |2002|) . Thus, 
Xd(LMC)=Xd(0) Zlmc=3.5x10-3 for Zlmc=0.5 Zq. We 
obtain an expression for the gas mass from the millimeter dust 
emission at 1.2 mm and for a dust temperature rd=20 K: 



M, 



Mr, 



158 



dust 



mJy 



(9) 



where ^'^'^^^^^^ is the cold dust millimeter continuum emission 
in mJy at 1.2 mm. The multiplicative factor varies from 236 for 
rd= 15 K to 119 for rd=25 K. This variation in temperature gives 
a maximum uncertainty on the millimeter masses of a factor of 
1.5. 

The resulting millimeter masses are listed in Table |4] 
All clouds but one, #12, are within 2cr of the median value 
Mvh/Mmm = 1.9 ± 0.4, where the en^or is the median absolute 
deviation. Cloud #12 is an outlier; it is compact, barely resolved 
by the observations and corresponds to a peak in the millimeter 
emission. We discard five molecular clouds because their esti- 
mated dust emission, after subtraction of the free-free emission, 
has a SNR lower than 2cr. These clouds are in the north-east of 
Nil, where the free-free emission is high (see Figure |9}. Their 
Mvh/Mmm ratio does not agree with the median value. 



5.4. Comparing the three mass estimates 

The observations provide three independent estimates of the 
masses of the molecular clouds in the Nl 1 region. Several pre- 
vious studies have compared the masses derived from CO lu- 
minosity with the virial masses for LMC molecular clouds. For 
the first time, in this paper, we add a third mass estimate based 
on the millimeter dust emission. Each mass estimate has a large 
systematic uncertainty, which we can constrain by quantifying 
what is required for the different methods to yield identical gas 
masses. 

The main uncertainty on the CO masses comes from the 
Xco-factor. The Xco-factor is the ratio between the true molecu- 
lar cloud mass and the CO luminosity. 



Xco - 



M,r 



i'CO(l-O) 



(10) 



As we pointed out in the introduction, in low metallicity 
environments the CO molecule is dissociated by the UV radi- 
ation from the nearby massive stars. Thus, the LMC Xco-factor 
must be larger than that measured in the Galaxy. However, in 
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Cloud 


^ 1.2mm 


5ff 


Sco 


■-^dust 


^mm 


Mvu- 


Mvi,-/M,„„ 


















mjy 






10^' 


Me 




1 


74±21 


13±6 


23+5 


27+23 


- 


125+11 


— 


3 


46±13 


8±4 


12+6 


23+15 


- 


31+18 


- 


4 


166±47 


21±10 


19+6 


123+49 


20+8 


28+11 


1.4+0.8 


5 


106±30 


36±17 


14+3 


53+35 


- 


72+15 


- 


7 


24±7 


0.2±0.1 


4+1 


19+7 


3+1 


8+3 


2.5+1.3 


8 


212±60 


36±17 


15+3 


159+62 


25+10 


42+14 


1.7+0.9 


9 


144±41 


64±30 


7+2 


68+51 


- 


49+20 


- 


10 


401±114 


125±59 


12+3 


261+123 


41+20 


50+8 


1.2+0.6 


11 


11±4 


2.3± 1.1 


3+1 


5+4 


- 


7+3 


- 


12 


142±40 


6+3 


6+2 


129+41 


20+6 


14+5 


0.7+0.3 


13 


324±92 


10±5 


22+5 


296+93 


47+15 


104+30 


2.2+0.9 


14 


284±81 


7±3 


40+8 


231+81 


37+13 


68+17 


1.9+0.7 


15 


93±27 


2.5+1.2 


5+2 


83+27 


13+4 


25+14 


1.9+1.3 


16 


54+16 


1.3+0.6 


5+1 


45+16 


7+3 


26+5 


3.6+1.6 


17 


24±8 


5.3+2.5 


5+1 


24+8 


4+1 


23+7 


6.0+2.9 


18 


84±26 


- 


22+5 


73+27 


12+4 


21+8 


1.8+1.1 


19 


38±12 


0.2+0.1 


7+2 


37+12 


6+2 


8+4 


1.4+1.0 



Table 4. Flux contribution to the total 1.2 mm emission for Nil molecular clouds. In the second row we list dust temperatures, 
obtained by a SED fitting (see section l53b . We separate this emission into 3 different contributions: free-free, CO(2-l) and cold 
dust emission. Millimeter masses are computed from the cold dust emission (see section l53] l. The error bar on dust masses (M,nn,) is 
derived from the error on 5 dust- It does not include any uncertainty on the dust temperatures. The last column shows the comparison 
between virial and millimeter masses. 



the LMC there is a wide range of values reported in the litera- 
ture, ranging froni a factor 2 to 6 higher tha n in th e Galaxy (e.g. 
Cohen etal."1988'; 'Israel" 1997t iFukui et al.ll200i iDobashi et all 
2008; Hughes et al. 2010). 

For the virial mass, we assume energy balance between the 
gravitational binding energy and turbulent kinetic energy. In do- 
ing this, we ignore external pressures as well as the magnetic 
energy. These contributions are difficult to quantify. The impor- 
tance of these energy terms can vary among the Nil clouds. 
Some of the molecular clouds could be influenced by the radia- 
tion pressure and winds from the near young massive stars. Other 
clouds could have a high external pressure because they are em- 
bedded within a bigger complex. Moreover, the energy terms are 
computed within the simplifying assumption of spherical geom- 
etry. To parametrize the difference between the true mass and 
the virial mass as estimated in section 15.21 we use the ffyir pa- 
rameter that is the ratio between the kinetic K and gravity W 
energy terms in the virial equations (Mc Kee & Zweibellll992h . 
For a molecular cloud with a density profile depending on the 
radius as p oc r"', the ffyir parameter is defined as 



2K 9 Rcrl 



\W\ 2GMa 






(11) 



The mass estimate based on the dust emission involves un- 
certainties on the dust temperature and the dust emissivity per 
hydrogen atom eH(v), which depends on the absorption coef- 
ficient of the dust Kiiiv) and the dust-to-gas mass ratio x^ (see 
section l53T l. We have computed the masses using a single dust 
temperature for each cloud but, on scales of a GMC, the dust 
has a range of temperatures. In Section ISTIJ we estimated the un- 
certainty on the dust mass associated with the choice of a single 
temperature to be less than 1 .5. The dust emissivity per hydrogen 
atom depends on two assumptions which are difficult to validate. 
We assume that the dust-to-gas mass ratio scales as the metal- 
licity and that the millimeter absorption coefficient /Cd(l-2mm) 
does not chan ge from the Milky Way to the LMC. This is a 
debated topic dPlanck Collaboration et al.ll201 id iGalliano et alJ 



1201 iL iDraine & Henslevll2012h and it is difficult to quantify un- 
certainties. To parametrize the uncertainties we introduce the 
factor TCdust = A'd(1.2mm)xci/[2.7xl0"^cm^g"*] where the refer- 
ence value is the value used to estimate the millimeter masses in 
section 5.3. This is related to the true molecular mass as follows. 



M,n,e = 



-C 1.2mm (dust) 



Mrr 



fil.2mm(rdust) /Cd(l-2mm)Xd TCdust' 



(12) 



where Xi.2mm(dust) is the dust luminosity at 1.2 mm per unit 
frequency, and M^^ is the millimeter mass computed in section 
5.3. 

The measured ratios between Lco/A^vii- and Mmm/Mvii- con- 
strain the three parameters introduced in formulae [TOJto [T2lif the 
three masses are to be identical. In the following equations we 
use the median values of these two ratios given in section lS!2l and 
15.31 and the error is the median absolute deviation of these val- 
ues. Note that these median values do not apply to all clouds. For 
a small number of clouds, the measured ratios are significantly 
diff'erent from the median values. 

First, the virial mass and the CO mass are equal if 



Xco = 8.8 ± 3.5 X 10^"cm"^ (K km s 



-')"«"> 



(13) 



and second, the virial mass and the millimeter mass are equal if 



•TCd, 



= 0.6 ± 0.2, 



(14) 



where the error bar does not include the uncertainty on the dust 
temperature. We assume that the mean temperature of 20 K ap- 
plies in average to the sample. 

From the previous equation, we estimate the product be- 
tween the millimeter absorption coefficient and the dust-to-gas 
mass ratio. 



A'd(1.2mm)xd = L5 : 



:0.5 xlO^^cm^ 



g ^v 



(15) 
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Fig. 10. The ratio between the product of the milUmeter ab- 
sorption coefficient Kd{l.2mm) and the dust-to-gas mass ratio 
Xd with the alpha parameter o-vir versus xj. This ratio is lin- 
early related to the ratio between the millimeter dust luminos- 
ity and the virial mass. The data represent the median values of 
the /(■(]( 1.2mm).5Cd/Q;vir ratio for molecular clouds located in three 
different environments with different dust-to-gas mass ratio, the 
Gould's Belt in the Milky Way, several regions in the SMC, and 
the Nil region in the LMC. The eiTor bar is the standard de- 
viation of each population. The values for the molecu lar clouds 
in the Gould's Belt and the SMC were extracted from lBot et al.l 
( |2007|) . The solid line is the Galactic measurement scaled to dif- 
ferent Xd values. The shaded area represents the uncertainty on 
the dust temperature (from 15 to 25 K). 



5.5. Comparison with Gould's Belt clouds in the Milky Way 

In Fig. [To] we compare our results for the N il clouds with 
earlier measurements done by iBot et al.l (l2007l) for clouds with 
comparable masses in the Gould's Belt (i.e. in the solar neigh- 
borhood) and the SMC. The figure shows the median ratio 
/('d(1.2mm)xd/Q'vir versus the dust-to-gas mass ratio Xd, for the 
molecular clouds in the three different environments. For a 
given dust temperature, the data points in Fig. [10] scale with 
the ratio between the dust luminosity at 1.2 mm and the virial 
masses (Eq. \T7\ . For the Galactic and SMC clouds, these val- 
ues were computed from the ratio bet ween the millime ter and 
virial masses listed in Tables 3 and 4 of'Bot e t aP (l2007h . Error 
bars are the standard deviation of the individual values. The solid 
line shows the Galactic value of /Cd(1.2mm)xd/a'vii scaled to dif- 
ferent gas-to-dust mass ratios. The pink region represents the 
uncertainty on the dust temperature, assuming a typical value 
of Td = 20 K which varies between 15 and 25 K. The me- 
dian value for the Nil clouds agree with the scaled Galactic 
value. This indicates that the observed difference in the mil- 
limeter dust luminosity to viiial mass ratio ratio between the 
Galactic and Nil clouds can be accounted for by the difference 
in the gas metallicities. This statement dos not apply to the SMC 
clou ds as reported by Rubio et al. (2004) and Bot et al. (2007, 
I2OIO ). The millimeter dust opacity k^, could be higher in the 
SMC than in the MiUcy Way. This interpretation is supported 
by the dif ference in the sub-mm/mm dust spectral energy dis - 
tributions ('Israel et al."2010^ iPlanck Collaboration et al.ll20TTcl) . 
[Oraine & Henslev (2012) have proposed that the enhanced mil- 
limeter dust emission observed in the SMC could be accounted 
for by dipolar magnetic emission from ferromagnetic particles. 



5.6. Discussion 

Our observational results constrain the product between the Xco- 
factor and the alpha virial parameter (Eq. [13), and the ratio be- 
tween the millimeter dust absorption coefficient and the virial 
parameter (Eq. \T5\ for molecular clouds in the Nil region. In 
this section, we discuss the most likely values for these three 
parameters. 

If the measured virial masses of the Nil molecular clouds 
are their true cloud masses, our Xco-factor is about 3 times 
larger than the canon ical value estimated for Galactic clouds 
(ISolomon et al.lll987l) . The effect of the undersampUng of the 
CO(2-l) observations will increase by a factor 1.5 the differ- 
ence between the averaged Xco -factor in the Milky Way and 
our estimate. Within uncertainties, it agrees with that found 
with the low -resolution (2!6 half-power be amwidth) NANTEN 
survey in the LMC by iFukui et aP (l2008l) . estimated from the 
virial and CO luminosity masses. But, it is about twice the value 
estimated by the MAGMA survey in the LMC (see Section ISTTT ) 
and that estimated by most of the studies of extragalactic GMCs, 
including GMCs from g alaxies in the local group (Bli tz et al. 
2007;'B olatto et al.l2008l) . as observed in Figure 7 in Lerov et al. 



(2011), where our Zco-factor estimate corresponds to cco = 
19 Mq pc-2 (Kkms-i)-'. The Xco-factor for the Nil clouds 
will agree with these previous estimates in the LMC if the avir 
parameter is a factor ~ 2 (see Eq.[T3]l. 

For molecular clouds in the LMC, I Wong et all (1201 ll) com- 
pared the ffvir paramete][f] with the CO luminous masses of the 
clouds (see the right panel of their Fig. 15). Using the canon- 
ical value of Xco, most of the LMC clouds have a^ii- > 1 
an d half of them gyjr > 2. Using the Xco-factor estimated 
by iHughes et al.l (1201 CJ) . half of the clouds have ayir > 1 and 
only a few gyjr > 2. For molecular clouds in the outer Galaxy, 
iHeveret aP (l200ll) found avh to be larger than 2 for clouds with 
CO masses smaller than 3 x lO"' Mq. Most of the Nl 1 clouds are 
more massive than this limit (see Table [2]i. How ever, in a recent 
study of ga lactic GMC s , defin e d previously by Solomon et al.l 
(11987!) and iHever et all (|2009|) . iDobbs et alj ([2011]) found that 
Galactic clouds, in the Nil mass range, have Oyt > 1 and about 
half of them have a cfvii > 2 (see their Fig. 1). This finding is 
in agreement with numerical simulations o f the turbulent ISM 
in the Milky Way presented in their paper. iDobbs et al.l (1201 ll) 
claim that there is no need for molecular clouds to be entirely 
bound to reproduce the observed star formation rate. Star forma- 
tion will occur in the densest part of the unbound clouds, where 
gravity wins over turbulence. 

For ffvir = 2, Equation[T5]yields a dust absorption coefficient 
for the Nil clouds of A-d - 0.86 cm g . This value is similar to 
the Galactic value presented in Section 1573] Therefore, we favor 
these values for the avir and k^ parameters because they do not 
require the dust emissivity, nor the dynamical state of the clouds, 
to differ between the LMC and the Milky Way. For avir = 2, we 
also get a Xco-factor akin to that estimated by the MAGMA sur- 
vey in the LMC. The ratio between the viiial parameter and the 
dust far-IR properties is plotted in Figure [10] The data does not 
allow us to separate the two parameters avir and k^. For Nl 1 with 
respect to the Galaxy we favor a variation of k^ associated with 
metallicity and no variation in avi,. However, this interpretation 
does not apply to SMC clouds. The figure shows that the ob- 
served /Cd(1.2mm)xd/avir values for the GMCs in the SMC are 
above than that estimated for the Galaxy and Nl 1, by a factor of 
3 to 4. 



This value depends on the assumed Xco-factor value. 
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6. Summary 

We have presented high sensitivity and spatial resolution obser- 
vations of the CO(2-l) line emission observed with SEST, for 
the Nil star-forming region in the LMC. Simultaneously with 
this CO transition, we observed the J=l-0 transition, which al- 
lowed us to measure the ratio between the CO(2-l) and CO(l-O) 
integrated emission for the entire map. This ratio has a median 
value of unity for the entire region. We identified 19 molecular 
clouds using the CPROPS algorithm which yields their CO lu- 
minosity, sizes, line-widths and virial masses. We also presented 
1 .2 mm continuum observations which we use to estimate mil- 
limeter dust fluxes. In this paper, and for the first time in the 
LMC, we compute molecular masses estimated from the mil- 
limeter dust luminosity (Xi.2mni(dust)) and compare them with 
the masses obtained from the CO luminosity (Leo) and virial 
theorem (Myir). The main results of this study are the following 
ones. 

- The correlations between CO luminosity, line-widths, sizes, 
and virial masses in the Nil clouds are in agreement with 
those found in earlier CO surveys of the LMC, sampling 
a larger set of clouds and a larger range of cloud masses 
(iFukui et al."2008t lWong et al.ll201 Ih . 

- The measured ratios Lco/^vir and Xi.2min(dust)/Mvir con- 
strain the Xco-factor and the dust emissivity per hydrogen 
atom at 1.2 mm as a function of the virial parameter ffyir 
(Mvii/Mtrue, the ratio between the kinetic and gravity energy 
terms in the viiial equations). 

- The comparison between the CO and virial masses 
yields a conversion factor of Xqo - 8.8 ± 3.5 x 
lO^^cm"^ (Kkm s"')"' c'^.'. The comparison between the 
virial and millimeter dust luminosity yields a product be- 
tween the gas-to-dust mass ratio x^ and the millimeter ab- 
sorption coefficient per unit of dust mass A:d( 1.2mm) of 
1.5 ±0.5 X lO-^cm^g-ia™-. 

- The comparison with previous studies of molecular clouds in 
the LMC suggests that Nil clouds are unbound, as observed 
in several Galactic molecular clouds within the same range 
of masses (iDobbset alJ 12011 ). The median values for the 

9 -1 

ttvir and Kd parameters are 2 and 0.86 cm g , respectively. 
These values do not require peculiar dust properties nor dy- 
namical clouds properties different from those observed in 
the Milky Way. 

- We do not find in Nl 1 a large discrepancy between the dust 
millimeter and virial masses, as reported for GMCs in the 
SMC by Rubio et al.1 (|2004) and Bot et alJ(l2007ll2010h . The 
ratio A:d(1.2mm)xd/avir in Nil is half of that measured for 
Gould's Belt molecular clouds, which can be accounted for 
by a factor two lower gas-to-dust mass ratio. This difference 
is the same as that observed between the gas metallicities. If 
the two samples have similar ayii median values, their dust 
far-IR properties should be also similar. 
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